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Summary. In this work we compare two semidiscrete schemes for the solution of 
hyperbolic conservation laws, namely the relaxation [JX95] and the Kurganov Tad- 
mor central scheme [KTOO]. We are particularly interested in their behavior under 
small time steps, in view of future applications to convection diffusion problems. 
The schemes are tested on two benchmark problems, with one space variable. 

1 Motivation 

We are interested in the solution of systems of equations of the form: 

Ut + fx(u) = Dp xx (u), (1) 

where f(u) is hyperbolic, i.e. the Jacobian of / is provided with real eigen- 
values and a basis of eigenvectors for each u, while p{u) is a non decreasing 
Lipshitz continuous function, with Lipshitz constant \i and D > 0. 

We continue the study of convection diffusion equations with the aid of 
high order relaxation schemes started in [CNPS06] for the case of the purely 
parabolic problem. 

In many applications, such as multiphase flows in porous media, p(u) is 
non linear and possibly degenerate. In these conditions, an implicit solution 
of the diffusion term can be computationally very expensive: in fact it may 
be necessary to solve large non linear algebraic systems of equations which, 
moreover, can be singular at degenerate points, i.e. where p(u) = 0. For this 
reason, it is of interest to consider the explicit solution of (1). This in turn 
poses one more difficulty. An explicit solution of (1) requires a parabolic CFL 
condition, that is stability will restrict the possible choice of the time step 
At to At < C(Ax) 2 , where Ax is the grid spacing. In other words, it may 
be necessary to choose very small time steps. But conventional solvers for 
convective operators typically work at their best for time steps close to a 
convective CFL, i.e. At < C Ax. When the time step is much smaller, they 
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exhibit a very large artificial diffusion of the form 0((Ax) 2r /At), where r is 
the accuracy of the scheme, see for instance [KTOO] . Clearly in these conditions 
artificial diffusion becomes very large for Zii — > 0. 

As a first step to the numerical solution of problem (1), we concentrate 
on semidiscretc schemes for the solution of the convective part of (1). Such 
schemes enjoy an artificial diffusion which depends weakly on At, and are 
therefore particularly suited for the solution of convection-diffusion equations. 

We will compare two semidiscrete methods for the integration of systems 
of hyperbolic equations. We are interested in the representation of solutions 
which can be characterized by strong gradients, and, in the degenerate case, 
even by discontinuities. Moreover, we are interested in comparing the behavior 
of the schemes for small values of At, and for such small values of the time step, 
we will investigate the resolution of discontinuous solutions and the behavior 
of the error in a few test problems. 

The schemes analyzed in this work are the Kurganov Tadmor central 
scheme proposed in [KTOO], and the relaxation scheme proposed in [JX95]. 
These methods discretize the equations starting from very different ideas, 
however they share some interesting properties. First of all, they are both 
semidiscrete schemes. Therefore they require separate discretizations in space 
and time, which is the key to the fact that artificial diffusion depends mainly 
on space discretization. Secondly, they are both Riemann solver free methods. 
The Kurganov- Tadmor scheme is based on a central approach: the solution of 
the Riemann problem is computed on a staggered cell, before being averaged 
back on the standard grid. In this fashion, the numerical solution is updated 
on the edges of the staggered grid, where it is smooth, and can be computed 
via a Taylor expansion, with no need to solve the actual Riemann problem. 
The relaxation scheme instead moves the non linearities of the convective term 
to a stiff source term, and the transport part of the system becomes linear, 
with a fixed and well known characteristic structure. Thus again there is no 
need to use approximate or exact Riemann solvers. 

For these reasons both schemes can be applied as black-box methods to a 
fairly general class of balance laws. 



2 Results 

For the Kurganov- Tadmor (KT) scheme we have followed the componentwise 
implementation of the method described in [KTOO] . The scheme is written in 
conservation form, with numerical flux: 



Fj+iMt) = \ [/(ut +1/2 (i)) + /(u7 +1/2 (i)) 



-a 3 +i/ 2 (t) (u+ +1/2 (t) 



J + l/2 



(*) 



(2) 



where 



(t) and Uj +1 / 2 (i) are the boundary extrapolated data, computed 
at the edges of each cell with a piecewise linear reconstruction at time t, and 
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a j+i/2{t) is a measure of the maximum propagation speed at the cell edge. 
This value for the case of systems of equations, in particular in the non convex 
case, must be carefully tuned, and it is the same for all components, when the 
scheme is implemented componentwise. 

On the other hand, the relaxation scheme requires an accurate choice of 
the subcharacteristic velocities A 2 . The relaxation system is 

du dv 

dt dx 

(3) 

dv . 2 du 1 

For a scalar conservation law, we take A 2 = max(|/'(u)|) as in [JX95], while 
for the Euler system of gas-dynamics we take A 2 to be the diagonal matrix 
with entries maxj (\uj—Cj\), max.,- ( | Uj | ) , maxj ( | Uj + Cj | ) . Here u is the velocity 
and c the speed of sound. We update these quantities at each time step, so 
that A 2 can be chosen as small as possible (in the paper [JX95] A 2 was chosen 
as a constant diagonal matrix but this results in a larger numerical diffusion) . 

Due to the diagonal form of A 2 , the convective operator is block diagonal 
with 2x2 blocks. Each block is independently diagonalized and we compute 
the numerical fluxes using a second order ENO reconstruction [HEOC87]. 

As e — > 0, the system (3) formally relaxes to the original conservation laws, 
provided the subcharacteristic condition holds, namely that {A 2 — (f'(u)) 2 ) is 
positive-definite. 

We use the second order Heun Runge-Kutta method for the time integra- 
tion of both the KT and the relaxation schemes. 
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Parabolic CFL 
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KT 
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2.03E-1 
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1.02E-1 
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7.58E-2 


7.66E-2 


2.04E-1 
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2.71E-2 
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8.22E-3 
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2.67E-2 


3.82E-3 


320 


2.29E-3 


2.29E-3 


7.62E-3 


1.03E-3 


640 


6.11E-4 


6.12E-4 


2.06E-3 


2.77E-4 


1280 


1.61E-4 


1.61E-4 







Table 1. Linear advection of a sine function. Errors in L 1 at t = 1. 



Table 1 shows the errors in the L 1 norm for the linear advection equation 
Ut + u x = with initial data u(x, 0) = sin(27rx). We use the standard con- 
vective CFL condition At = CAx and the parabolic CFL, At = C(Ax) 2 . Wc 
note that the errors are almost the same for the two schemes for the convective 
CFL, while the relaxation scheme seems superior for the parabolic CFL. 
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A key requirement for a numerical scheme for conservation laws is the 
ability to pick the entropy solution in non-convex problems. Here we show 
a Riemann problem for the non-convex flux f(u) = (u 2 — l)(u 2 — 4)/4, as 
in [KTOO]. The Riemann problem breaks into two shocks connected by a 
rarefaction wave. The results are shown in Figure 1. Clearly both schemes are 
able to resolve the correct discontinuities and they have approximately the 
same resolution, the KT scheme being slightly less diffusive. 




Fig. 1. Non-convex flux. Kurganov-Tadmor and relaxation schemes, with n — 100 
(left) and n = 200 (right). 

Figure 2 shows the density component of the Lax Riemann problem in 
gas dynamics. The accurate choice suggested above for the matrix A 2 in the 
relaxation system yields a slightly higher resolution than KT. 




Fig. 2. Lax Riemann problem, density component. KT (dashed) and relaxation 
(solid line) with n = 100 (left) and n = 200 (right), where a detail of the density 
peak is shown. 
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Concluding remarks 

We have compared two semidiscrete schemes for conservation laws. We find 
that although the schemes are constructed with very different philosophies, 
they yield comparable results on some significant test problems. We think 
that the relaxation scheme is slightly more robust, since it results from the 
relaxation of a viscous profile, provided the subcharacteristic condition is sat- 
isfied. Also, the actual errors obtained with a parabolic CFL in Table 1 seem 
to favour the relaxation scheme. 

We also wish to mention higher order extensions of the schemes stud- 
ied in this work: namely the third order central upwind scheme described 
in [KNP01], endowed with a more carefully crafted artificial diffusion with 
respect to [KTOO] and the third order extension of the relaxation scheme pro- 
posed in [Sea06]. 
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